VT symmetry and spontaneous symmetry breaking in a microwave billiard 
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We demonstrate the presence of parity-time (VT) symmetry for the non-Hermitian two-state 
Hamiltonian of a dissipative microwave billiard in the vicinity of an exceptional point (EP). The 
shape of the billiard depends on two parameters. The Hamiltonian is determined from the measured 
resonance spectrum on a fine grid in the parameter plane. After applying a purely imaginary 
diagonal shift to the Hamiltonian, its eigenvalues are either real or complex conjugate on a curve, 
which passes through the EP. An appropriate basis choice reveals its VT symmetry. Spontaneous 
symmetry breaking occurs at the EP. 
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Introduction. Flat microwave cavities, so called "mi- 
crowave billiards" are analogues of quantum billiards [1]. 
They are an experimental test ground for the properties 
of the eigenvalues and eigenfunctions of quantum sys- 
tems [2]. In the present Letter we demonstrate theoret- 
ically and experimentally that they can also be used to 
study dissipative quantum systems which have a parity- 
time (VT) symmetry, that is, are invariant under the 
simultaneous action of a parity operator (P) and a time 
reversal operator (T) after a suitable width-offset. 

Strong interest in PT-symmetric quantum systems 
was initiated in 1998 by Ref. [3] demonstrating that a 
non-Hermitian Hamiltonian H has real eigenvalues pro- 
vided it respects VT symmetry, i.e. [PT,7-L] = and 
has eigenvectors that are also VT symmetric. In [4] this 
statement was generalized to systems invariant under an 
antilinear operator T' = UT, where U is unitary. The 
observations in Refs. [3|, |4[ led to a reconsideration of 
the necessity of the Hermiticity axiom for quantum ob- 
servables [5| . Depending on external parameters the VT 
symmetry of the eigenvectors may be spontaneously bro- 
ken, i.e. they cease to be eigenvectors of PT ', although H 
still commutes with PT [3|, |5[. As a result, the eigenval- 
ues of H are no longer real, but rather become complex 
conjugate pairs. This phase transition occurs at an ex- 
ceptional point (EP) [6] defined as the coalescence of at 
least two eigenvalues and the corresponding eigenvectors. 

Close to an EP, non-Hermitian but T^T-symmetric 
Hamiltonians are capable to speed-up quantum evolution 
processes [7j]. During the last years optical wave-guide 
systems with fine-tuned (PT-symmetrically balanced) 
gain and loss regions, i.e. with active (laser pumped) 
and absorptive components, have been investigated the- 
oretically [8[. This led to experimental setups which 
consist of one wave-guide with absorption and another 
one which is either optically pumped (active VT sym- 
metry) [9| or lossless (passive VT symmetry) [10|. In 
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FIG. 1: Two typical spectra |£ii(/)| 2 for B = 38 mT and 
the parameter settings (s,S) = (1.66,41.79) mm (solid line), 
and (s, S) = (1.99,41.76) mm (dashed line), respectively, 
in the vicinity of the EP, which is located at (sep,#ep) = 
(1.72 ± 0.01,41.78 ± 0.01) mm. Inset: top view (to scale) of 
the microwave billiard. The setup is described in the text. 



optical systems the VT phase transition has been stud- 
ied theoretically [11] as well as experimentally [9|, llQj . 
Recently it was observed in a non-Hermitian system 
with active VT symmetry consisting of one amplifying 
and one attenuating LRC circuit [12|. Further theo- 
retical studies of VT symmetry based effects concern 
spectral singularities [13], lasers at threshold [14J, [l5j, 
coherent perfect absorbers [15[, unidirectional invisibil- 
ity induced by T^T-symmetric periodic structures [16[ as 
well as Bloch oscillations in T^T-symmetric lattice struc- 
tures [17] and optical tachyons in PT- symmetric optical 
graphene-like structures [18j]. The ^-matrix formalism 
for PT-symmetric systems was analyzed recently [15lll9j. 
Experiment. We will show that configurations with 
passive VT symmetry, including a VT phase transition, 
are observable in a dissipative microwave billiard, whose 
shape depends on two parameters. The experiments were 
performed in the vicinity of an EP in the parameter 
plane. The setup was similar to that used in [20|. There 
the first experimental evidence of EPs and of geometric 
phases was provided. Recently, these experiments have 



been extended to systems with induced T violation [21 1. 
The obtained data are used for the present work. The 
experimental setup (see the inset of Fig. [1]), described 
in [2l|, consisted of a flat, 5 mm high microwave res- 
onator of circular shape with 125 mm radius. It was 
divided into two approximately equal parts by a 10 mm 
thick copper bar, with an opening of 80 mm. The cou- 
pling between the electric field modes in each part via the 
opening was varied with a copper gate that had a tilted 
bottom. It was inserted into the resonator through a 
slit in its top and could be moved up and down. The 
bottom plate had a notch to enable a complete clos- 
ing of the gate. The lifting s defines one parameter, 
where (no coupling) < s < 9 mm (maximal coupling). 
The second parameter is provided by the displacement S 
(with respect to the cavity center) of a 5 mm high semi- 
circular Teflon piece of radius 30 mm inserted into the 
left part of the cavity. The parameter plane (s, 5) was 
scanned with the help of two micrometer stepper mo- 
tors, which moved the gate and the Teflon piece in steps 
of As = A5 = 0.01 mm. Furthermore, a 5 mm high 
cylindrical ferrite of radius 2 mm denoted by F in Fig. [1] 
was inserted into the right part of the cavity. To induce 
T violation the ferrite was magnetized with an external 
magnetic field B of strength < B < 90 mT, applied per- 
pendicularly to the billiard plane [2l|, [22| . Two pointlike 
wire antennas 1 and 2 reached into the cavity, one into 
its left, the other into its right part. A vectorial network 
analyzer Agilent PNA 5230A emitted microwave power 
into the resonator via one antenna a and received an 
output signal either at the same or at the other antenna 
b. It determined the amplitude and phase of the output 
signal relative to the input signal, thus yielding the four 
complex elements Sba(f), { a ,b} G {1,2} of the scatter- 
ing matrix S(f). They were measured for each setting of 
the parameters (s, S) as a function of the excitation fre- 
quency / in steps of A/ = 10 kHz. The frequency range 
of 40 MHz was determined by the spread of the reso- 
nance doublet under consideration. In that range the 
electric field vector is perpendicular to the top and bot- 
tom plates. Therefore, the Helmholtz equation is math- 
ematically identical to the Schrodinger equation of the 
quantum billiard [l|, [2[. Thus, the results of this Letter 
also apply to the associated quantum system. Figure Q] 
shows two reflection spectra measured for B = 38 mT, 
one for s < Sep (solid line) and one for s > Sep (dashed 
line). For the former we observe a close encounter of 
the resonance positions, for the latter one of the widths. 
In fact, we observe this feature of the resonance spectra 
along a curve in the parameter space which crosses the 
EP independently of the choice of the external magnetic 
field. Exactly at the EP the resonance shape is that of a 
second order pole in addition to the first order poles [29]. 
As demonstrated in the following, along this curve the 
effective Hamiltonian has the form of a PT-symmetric 
Hamiltonian after a suitable basis transformation. 



Neighboring resonances are situated about 250 MHz 
away from the doublet under consideration. Con- 
sequently, the effective Hamiltonian H e ^ is two- 
dimensional. It is determined for every setting of s,5,B 
by fitting an analytic expression for the S- matrix 22|, [23| 
originally derived in the context of nuclear reaction the- 
ory and extended to microwave resonators in [25| , 
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to the measured one. The effective Hamiltonian 
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is obtained by evaluating the integrals entering 
Eq. (4.2.20b) of Ref. [H| and Eq. (4) of Q. The matrix 
1L is the unit matrix and H is the Hermitian two-state 
Hamiltonian of the closed resonator or, equivalent ly, the 
quantum billiard. It includes the coupling of the ferro- 
magnetic resonance to the rf magnetic field in case of 
a nonzero external magnetic field B and thus takes ac- 
count of the violation of T invar iance. The matrix el- 
ements Wafj,, W\)^ couple the resonator modes /i = 1,2 
to the antenna states {a, b} G {1,2} and the sum over 
c includes the fictitious channels, which describe dissipa- 
tion in the walls of the resonator and the ferrite [24J, |26| . 
These matrix elements are real and frequency indepen- 
dent in the considered range. Thus, for a vanishing ex- 
ternal magnetic field H e ^ is given by Eq. ((2]) with real 
matrix elements W a ^ and a time reversal invariant H, 
[T, H] = 0. Here, T is the antilinear operator of complex 
conjugation, whence H is real symmetric. As a conse- 
quence, the scattering process is also time-reversal in- 
variant, TST = S^ [27j, and the S matrix is symmetric, 
S a b = Sba- T violation is induced with a nonvanish- 
ing 5/0. In this case, the coupling of the resonator 
states to the ferromagnetic resonance is complex. Conse- 
quently, the evaluation of the integrals entering Eq. (4) 
of Ref. [22j yields a complex and antisymmetric contri- 
bution, leading to a complex Hermitian Hamiltonian H 
in Eq. (J2j), which is not invariant with respect to T. Lack 
of reciprocity, i.e. S a b 7^ Sb a , in the measured spectra is 
the signature for T violation [2l|, [22| . In the T- invariant 
case i^ eff is complex symmetric, otherwise it is nonsym- 
metric. To simplify notation we express it most gener- 
ally in terms of the unit matrix 1L and the Pauli matrices 
<? — fox: cFy> cfz) as ( see Eq. (1) of Ref. [21(), 
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(3) 
All entries of H eS are complex, and depend on (5, (5), 
but not on / in the considered frequency range. Fitting 
the S matrix Eq. ([T]) to the measured one determines 
H ef[ and W\^ W^/z U P to common real orthogonal basis 



transformations. As in [21( we choose the basis such that 
the ratio of the off-diagonal elements of H eE equals 
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with a real T- violation parameter r. Maximal T vio- 
lation occurs for r = ±7r/4. In the T-invariant case 
H^2 = [22] one has r = 0. The transformation to this 
basis is achieved with Oq = e l ® 0(Jy where, in terms of the 
components of h, tan(2<l>o) = S (/13//12) /^ (^1/^2)- The 
fitting procedure was tested thoroughly in [21], via com- 
parison with the method used in [20] to determine the 
complex eigenvalues of H eS . With Eq. (J3j) the complex 
eigenvalues Ej = fj — iTj/2 of H eE are given as 
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At an EP the radicand V vanishes, so that 

with |!R/i| 2 , \^sh\ 2 nonvanishing, and Wi • ^sh 

EP was found by proceeding as described in [2l|. It was 

confirmed by determining the geometric phases gathered 

by the eigenvectors of H e ^ on encircling the EP. 

The real parts of the eigenvalues yield the positions, 
the imaginary parts the widths of the resonances in the 
measured spectra \S a b (/)| 2 - Since the system is dissi- 
pative, the widths Ti^ > are nonvanishing and, con- 
sequently, the eigenvalues are complex in the whole pa- 
rameter plane. Yet, we will demonstrate in the following, 
that there exists a curve in the parameter plane along 
which the eigenvalues £1^ = ^1,2 + ifii + ^2)/^ of the 
non-Hermitian Hamiltonian Ti = H e ^ + z!L(ri + IV) /4 
are either real or complex conjugate. This entails [3|45| 
that Ti has a VT symmetry along this curve, which is of 
so called passive type. The situation is similar to that 
in the experiments with optical wave-guides without ac- 
tive laser pumping [10j and in Bose-Einstein condensates 
with leakage and without injection [28|- A pair of real 
eigenvalues typical for T^T-symmetric systems would cor- 
respond to resonances of zero width [13| . It is obtained 
only after a suitably chosen width-offset. This can be 
avoided by a pumping mechanism [9|, LL2|, [l9j. However, 
the eigenvectors and thus the physics of VT symmetry 
and of the phase transition to spontaneously broken VT 
symmetry at an EP is not affected by this shift. Fur- 
thermore, the dissipation and, consequently, the shift de- 
pends only marginally on the parameters s and <5, and 
thus the excitation frequency. Accordingly, we will con- 
sider Ti instead of i7 eff . Note that T does not commute 
with the non-Hermitian Ti, not even when the internal 
Hamiltonian H is T invariant, as is the case for B = 0. 
Yet there is a curve in the parameter plane on which Ti 
is VT symmetric. 

Experimental results. We located the curve in ques- 
tion by plotting for a given B the difference E\ — E2 = 



£ 1 — 82 = 2yV in the parameter plane. The left and right 
panels of Fig. [2] show the distances of the eigenvalues for 
B = and B = 38 mT, respectively. The darker the color 




1.6 1.8 

s (mm) 



FIG. 2: (Color online) Differences of the complex eigenvalues 
-Ei ,2 = /i,2 — iTi^/2 in an area of the parameter plane (s, 6) 
around the EP located at (s E p,£ep) = (1.68 ± 0.01,41.19 ± 
0.01) mm for B = (left panel) and at (s E p,fcp) = (1.72 ± 
0.01,41.78 ± 0.01) mm for B = 38 mT (right panel). The 
darker the color the smaller is the respective difference. 

the smaller is the respective distance. The color scale was 
chosen such that only distances |/i — /2I < 3 MHz and 
|Ti — T2 1 < 0.35 MHz are shown colored. We observe, 
that in both panels the distances of the real parts of the 
eigenvalues are small to the left side of the EP, i.e. for 
s < Sep, those of the imaginary parts for s > Sep and 
both vanish at the EP, (s,5) = (sep,#ep)- Along the 
curve of darkest color they are vanishingly small. Fig- 
ure [3] shows the three parts of V defined in Eq. (J5]), 
|3?/i| 2 , \^sh\ 2 and $lh • S/i, separately along this curve, 
in the left panel for B = and in the right one for 
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FIG. 3: Radicand V defined in Eq. JS} for the parameter 
values (s, 6) along the curve of darkest color in Fig. [2] Shown 
are \^fth\ 2 (crosses), \^sh\ 2 (circles) and ?fth • Ssh (triangles), in 
the left panel for 5 = 0, in the right one for B — 38 mT. To 
obtain a dimensionless h, it was divided by \h±\. 

B = 38 mT. Note the continuity of the curves, which 
demonstrates the precision of the measurements and of 
the determination of H eE . In fact, each point was ob- 
tained from an independent fitting. For both values of 
5, Wi • ^sh (triangles) is vanishingly small, that is the 
radicand V is approximately real along the dark curves 
in Fig. [21 The curves \Wi\ 2 (crosses) and \^sh\ 2 (circles) 
cross at the EP. For s > Sep the former is larger than 
the latter; consequently \[T> and thus the eigenvalues of 
Ti are real. They change into a pair of complex conju- 
gate eigenvalues for s < Sep- The transition takes place 
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FIG. 4: The parameter r in units of maximum T violation 
7r/4 along the curve of darkest color in Fig. [2] for B — 38 mT. 



at the EP. This behavior is reminiscent of that of a VT- 
symmetric Hamiltonian complemented by a spontaneous 
breaking of the VT symmetry of its eigenvectors at the 
EP. 

Interpretation of the experiment. To reveal the VT 
symmetry of % along the curve $lh ■ ^h = we search for 
the basis transformation which yields the parity operator 
as P = cr x . Most generally, in this basis a Hamiltonian 
that commutes with PT can be written in the form 
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The eigenvalues of H PT are real when its eigenvectors 
also commute with PT ', complex conjugate otherwise. If 
T invariance is violated Hfy in Eq. (|3j) and r in Eq. ((4]) 
are nonvanishing. The T-violation parameter r varies 
with the parameters s and £, because the wave func- 
tions [20, [30| and thus the electromagnetic field at the 
position of the ferrite change with the opening size s 
of the slit and the position 5 of the Teflon piece. In 
Fig. |4] we show for B = 38 mT its change along the 
curve characterized by $lh ■ ^sh = 0. On this curve H can 
be brought into the form of Eq. ((6]). The transforma- 

tion U = e- W>*. = I ex P("W2) \ tumg 

y exp(zr/2) J 

% into a complex symmetric matrix. For a T-invariant 
scattering process r = and 17 = 11, because then % is al- 
ready complex symmetric [see Eq. ([3])]. If and only if V is 
real, i.e. $lh-^sh = 0, there is a basis transformation that 
brings Ul-LU^ to the form of Eq. ([6]). This is achieved with 
the real orthogonal transformation O = e l ® (Jy , where 
tan(2<£) = 9/ii/ [S7i3 cos r]. The transformation yields 
Eq. © with A = (/i + / 2 )/2, B = Q?^/ [sin (2$) cost], 
C = sRfti/ [cos (2$) cos r] and D = 0. 

Conclusions. We experimentally identified a curve in 
the parameter plane along which the eigenvalues of % are 
either real or complex conjugate. There, the radicand V 
is real, i.e. Wi • Ssh = 0. This relation is independent of 
the basis representation of H. We have specified the ba- 
sis transformation, which brings T-L to the form of Eq. (J6j) 
along the curve. In that basis, VT symmetry corresponds 
to invariance with respect to PT. If the scattering sys- 
tem is T invariant, that is for B = 0, the basis change 
is achieved with the real orthogonal transformation O, 



i.e. [PT, OHO T ] = 0. In distinction to previous exper- 
iments [Ol, ll0|, LL2| , the present data include T violation 
of the scattering process, induced with an external mag- 
netic field 5/0. In that case, the unitary basis trans- 
formation OU brings H to the form of Eq. (|6j), that is, 
[Pf,6imU^6 T ) =0.^ Thus, H is invariant under the 
antilinear operator U'T, where U' = ffi T POU^ is uni- 
tary, in accordance with [4]. As predicted, the change 
from real eigenvalues for s > s EP to complex conjugate 
ones for s < s EP is accompanied by a sponteneous break- 
ing of VT symmetry of the eigenvectors of OUl~LWO T at 
the EP, that is, they cease to be eigenvectors of PT [3|,[5|- 
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